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IMAGE PROCESSING METHOD FOR AUTOMA TI C ADAP TATION OF 3-D DEFORMABL E MODEL ONTO 
A SUBSTANTIALLY TUBULAR SURFACE OFTL 3-D OJETfiCT 


Description 
5 Field of the Invention 

The invention relates to an image processing method for automatic adaptation 
of a 3-D deformable mesh model onto a substantially tubular surface of an object in a 
medical image. The invention applies to methods of segmentation of a three 
dimensional tubular object in a three dimensional image, which often comprises an 

10 operation of fitting a three-dimensional discrete deformable model onto said three- 
dimensional tubular object. The invention further relates to medical imaging apparatus 
or systems and to program products for processing medical three-dimensional images 
produced by those apparatus or systems, for the segmentation of objects that are body 
organs such as colon or arteries in order to study or detect organ pathologies. The 

15 invention finds a particular application in the field of medical imaging methods, 
program products and apparatus or systems. 

Background of the Invention 
A technique of modelization of a 3-D object is already disclosed by H. 
DELINGETTE in the publication entitled "Simplex Meshes: a General Representation 

20 for 3D shape Reconstruction" in the "processing of the International Conference on 

Computer Vision and Pattern Recognition (CVPR'94), 20-24 June 1994, Seattle, USA". 
In this paper, a physically based approach for recovering three-dimensional objects is 
presented. This approach is based on the geometry of "Simplex Meshes". Elastic 
behavior of the meshes is modeled by local stabilizing functions controlling the mean 

25 curvature through the simplex angle extracted at each vertex (node of the mesh). Those 
functions are viewpoint-invariant, intrinsic and scale-sensitive. A Simplex Mesh has 
constant vertex connectivity. For representing 3-D surfaces, Simplex Meshes, which are 
called two-Simplex Meshes, where each vertex is connected to three neighboring 
vertices, are used. The structure of a Simplex Mesh is dual to the structure of a 

30 triangulation as illustrated by the FIG.l of the cited publication. The contour on a 

Simplex Mesh is defined as a closed polygonal chain consisting of neighboring vertices 
on the Simplex Mesh. Four independent transformations are defined for achieving the 
whole range of possible mesh transformations. They consist in inserting or deleting 


WO 2004/063988 


2 


PCT/IB2004/000063 


edges in a face. The description of the Simplex Mesh also comprises the definition of a 
Simplex Angle that generalized the angle used in planar geometry; and the definition of 
metric parameters, which describe how the vertex is located with respect to its three 
neighbors. Dynamic of each vertex is given by a Newtonian law of motion. The 
5 deformation implies a force that constrains the shape to be smooth and a force that 

constrains the mesh to be close to the 3-D object. Internal forces determine the response 
of a physically based model to external constraints. The internal forces are expressed so 
that they be intrinsic viewpoint invariant and scale dependant. Similar types of 
constraints hold for contours. Hence, the cited publication provides a simple model for 

10 representing a given 3-D object. It defines the forces to be applied in order to reshape 
and adjust the model onto the 3-D object of interest. 

Summary of the Invention 
In medical images, it is often required to segment a tubular organ like the colon 
or an artery. A segmentation based on discrete deformable models allows to extracting 

15 clinical parameters of the studied organ like the diameter or the volume. Problems arise 
when the discrete deformable model, whether of the kind called two-Simplex Mesh, 
triangular Mesh or of any other kind of active contour Models, must fit such a tubular 
organ that presents numerous bends. These bends may have different curvatures, 
defined by small or large curvature radii. 

20 It is an object of the present invention to propose an image processing method 

of creating a deformable model, with 2-simplex meshes or triangular meshes or any 
other kind of meshes, having a tubular structure for fitting a previously defined 3-D path 
composed of a set of ordered points. This 3-D path may correspond to the centerline of 
an object of interest, i. e. a tubular organ showing bends. Thus, said 3-D path may 

25 present all kinds of curvatures. It is a further object of the present invention to propose 
such an image processing method for automatic adaptation of the mapping of the 3-D 
deformable tubular model onto the 3-D surface of the tubular object of interest showing 
different bends. The object of interest may be represented in a gray level 3-D image. 
This method comprises: computing a 3-D path that corresponds to the centerline of a 

30 tubular object of interest to segment; defining segments on said 3-D path; creating an 
initial straight deformable cylindrical mesh model, of any kind of mesh, with a length 
defined along its longitudinal axis equal to the length of the 3-D path; dividing this 
initial mesh model into segments of length related to the different segments of the 3-D 
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path; and computing for each segment of the mesh, a rigid-body transformation that 
transforms the initial direction of the mesh into the direction of the related segment of 
the 3-D path, and applying this transformation to the vertices of the mesh corresponding 
to that segment. 

5 Now, a key problem lies in the potential apparition of self-intersections of the 

mesh model when, due to the bends of the 3-D path, the direction between two 

consecutive segments of the 3-D path changes quickly. 

It is an object of the present invention to propose such an image processing 

method capable of limiting self-intersections of the deformable mesh model. For 
10 limiting self-intersections, the proposed method does not use a unique transformation 

for each segment, but instead, the proposed method comprises computing 

transformations related to following segments, which transformations, are blended in 

between two consecutive segments. This technique is called 3-D rigid body 

transformation blending. When rotations are used for following the directions of 
1 5 segments, a linear interpolation can be used between two rotations for rotation blending. 

This linear interpolation, from one 3-D rigid body transformation to another one can be 

solved using the formalism of Quaternions. 

An advantage of this technique of 3-D rigid body transformation blending is 

that the creation of the mesh model according the proposed method does not need to 
20 impose a kind of mesh model rather than another kind. The proposed method can be 

applied to mesh model of any kind: 2-simplex mesh, triangular mesh or any other kind 

of transformable model. 

Linear interpolation of 3-D rigid transformation from one segment to the other 

does not always suffice to completely avoid self-intersections. Clearly, such self- 
25 intersections also depend on the relation between the local curvature of the 3-D path and 

the desired radius of the created discrete deformable model. If the latter is larger than 

the local radius of curvature, then self-intersections occur. 

It is an object of the present invention to propose such an image processing 

method capable of avoiding possible self-intersections in the regions of the tubular 
30 deformable model showing bends. The method of the present invention comprises 

modulating the radius of the cylindrical deformable model according to the local 

curvature. First, since the target 3-D path can be nothing more than a mere list of points, 

local curvature can be approximated. Then, the radius of the mesh model is restricted in 
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function of the curvature. However, in order to avoid sharp radius changes from one 
segment of the model to the other, a radius modulation technique is applied to the radius 
of the mesh model, such as linear blending or bi-cubic spline interpolation from one 
radius to the other. 

5 An advantage of this technique of radius modulation is that the surface of the 

mesh model, in the parts with restricted radius obtained by radius modulation, does not 
present any discontinuity. The surface passes smoothly from a part of a given radius to 
another part of restricted radius. Hence the mesh model does not appear to be deformed. 
Another advantage of this technique of radius modulation is that it permits an improved 

10 segmentation of the complex organ of interest in regions where numerous folding up 
may make visualization or following of all the parts of this organ difficult. The 
visualization is improved. 

Another advantage of this technique of radius modulation is that the creation of 
the mesh model according the proposed method does not need to apply truncation to the 

1 5 self-intersecting parts of the bends, since the radius of the mesh model is restricted. The 
avoided alternative of truncation is a kind of operation that those skilled in the art know 
to apply to the parts of the bends that self-intersects. Such truncation could be a solution 
for avoiding self-intersections, but presents drawbacks such as suppression of 
information in the regions of curvatures and un-natural abrupt folds, which the present 

20 invention has for advantage to avoid. 

Still another problem may be lack of continuity control of the above-cited 
transformation, which can be seen as mesh torsion in the model. 

It is a further object of the present invention to propose such an image 
processing method capable of minimizing Mesh Torsion, Mesh torsion is minimized 

25 when the distance between two consecutive rotations of the rigid-body transformations 
is minimal. The proposed method comprises computing the minimal 3-D rotation from 
the initial mesh direction to the target segment. This rotation is defined with an axis 
parameter and with a rotation angle parameter. The proposed method comprises 
computing these parameters iteratively from one segment to the other so that the new 

30 rotation for a current segment is computed as a composition of the found rotation for the 
previous segment and the minimal rotation from the previous and the current segment. 

The main features of the proposed image processing method are claimed in 
Claim 1 . The method of the invention can favorably be applied to the segmentation of 
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three-dimensional medical images of tubular body organs. The invention also relates to 
a medical examination imaging apparatus having 3-D image processing means. The 
medical imaging apparatus may be any 3-D medical examination apparatus having 
imaging means. The invention further relates to a program product or a program 
5 package for carrying out the method. 

Brief Description of the Drawings 
The invention is described hereafter in detail in reference to the following 
diagrammatic drawings, wherein: 

FIG.l illustrates the step of mesh bending segment by segment, based on a 
1 0 predetermined path of ordered points; 

FIG.2 illustrates the variation of the angle of the local tangent u(s) to the 3-D 
path withrespect to the z-axis of the referential of the initial tubular model; 

FIG.3 A and FIG.3 B illustrate respectively mesh creation without and with 
linear transformation blending, in circle views; 
1 5 FIG. 4A and FIG.4 B illustrate respectively mesh creation without and with 

linear transformation blending, in simplex mesh views; 

FIG.5 shows an example of bi-cubic spline interpolation on computed 
curvatures; 

FIG.6A shows an example of mesh model creation with linear transformation 
20 blendingandtofSon minimization, from the path of a centerline of a colon; 

FIG.6B shows the example of FIG.6A related to a mesh model following the 
path of a centerline of a colon, further modified by a step of radius modulation; 

FIG.7A and FI G.7Bi llustrate mesh creation using the minimal rotation 

^ . 

between the local tangent u(s) and the z-axis without and with incremental rotation R(s), 
25 in mesh views; 

F TO.8A shows an initial tubular mesh model; FIG.8B shows a synthetic 
tortuous example of mesh creation with linear transformation blending and with 
incremental rotation R(s) leading to torsion minimization; FIG.8C shows the same 
example as _FIG.8Bjn an other orientation; and FIG.8D shows the same example as 
30 FIG.6B with a finer resolution of the 3-D path; 

FIG.9 illustrates a medical viewing system having means for carrying out the 
image processing method. 

Description of the preferred Embodiments 
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The invention relates to an image processing method to be applied for example 
to a three-dimensional (3-D) digital image represented in gray levels. The image may 
represent the noisy three-dimensional surface of an organ called object of interest. In 
order to provide the user with a better view of the object of interest, for instance with 
5 respect to the background, this object is segmented. The segmented image permits the 
user of better studying or detecting abnormalities or diseases of the organ. In the present 
example, the object of interest is a tubular organ, for instance the colon or an artery or 
the bronchi. The way the three-dimensional image is acquired is not part of the 
invention. The segmentation method could be applied to three-dimensional digital 

1 0 images of organs that can be acquired by ultrasound systems or X-ray apparatus or by 
other systems known of those skilled in the art. 

After the acquisition of the three-dimensional image representing the three- 
dimensional tubular object of interest, said image is segmented. The image 
segmentation method is based on the utilization of 3-D deformable models, called active 

15 contours. According to the invention, any technique of creating a 3-D deformable model 
can be used without restriction. The segmentation operation consists in mapping the 3-D 
deformable Model onto a 3-D substantially tubular object of interest, such as a colon or 
artery. In this example, the pbject of interest shows a complex tubular shape comprising 
numerous bends. 

20 In the field of active contours, an initial mesh model has to be provided. Even 

if it is always possible to start from any arbitrary shape of the mesh model, it is more 
robust and faster to start with a mesh model whose shape is close to the desired shape of 
the organ to be segmented. According to the invention, creating an initial mesh model 
of the kind called simplex mesh, triangular mesh or any other kind of mesh model is 

25 proposed, with a tubular structure fitting a 3-D path formed of a set of ordered points. 
This method is based on the creation of a long straight cylinder, which is in turn bended 
or deformed to fit the 3-D path. 

The difficulty lies in the operation of deforming the straight initial tubular 
deformable Model appropriately in order to finally map correctly the entirety of the 

30 surface of the complex shape of this tubular body organ considering the complexity of 
said shape. 

1) Creating a tabular mesh model: 
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Referring to FIG.l, the segmentation of a tubular structure, like a blood vessel 
or like the colon, comprises to first create a 3-D path P. This 3-D path may 
advantageously be a centerline of the organ corresponding to said tubular structure. The 
creation of this 3-D path P may be performed using path-tracking techniques known to 
5 those skilled in the art for providing the set of ordered points related the tubular 

structure. According to the invention, the 3-D path P is converted into a tubular mesh 
model. The initial mesh model is a cylindrical mesh model. This initial cylindrical mesh 
model is deformed towards the tubular organ structure data. For this, functionality is 
needed in order to initialize a mesh model from such a 3-D path P, instead of initializing 
10 a mesh model directly with an object surface as in the prior art publication [Delingette]. 

As mentioned in the introduction, any application aiming at segmenting tubular 
structures might benefit from an initial mesh model having a tubular shape. Since path- 
tracking tools are already known to those skilled in the art, they may be used to 
determine a line of points related to the object of interest such as the centerline of the 
15 tubular object of interest to be segmented. Such a 3-D path P, or centerline, can guide 
trustfully the initialization process of the mesh model, leading not only to a faster but 
also to a more robust segmentation tool for tubular structures. 

According to the invention, the problem to be solved is the creation of the 
tubular mesh model fitting the tubular organ to be segmented. The different inputs of 
20 the problem are: 

1) a sorted list of points lying along the 3-D path P. No assumptions are 
required yet on the regularity and the spacing of these points, but it will be 
demonstrated further that such constraints can really help in obtaining a smooth mesh 
model. 

25 2) the radius r of the cylinder, and 

3) the resolution of the cells, i.e. the "number of cell subdivisions" . 
The natural output is a mesh structure. 

Referring to FIG.l, a functionality for creating the cylinder basic form is 
proposed. This functionality consists in creating along the z-axis of a predefined 
30 referential Ox, Oy, Oz, a set of points lying on circular sections of the initial cylindrical 
mesh model, then linking the sets of points all together to create the simplex mesh 
structure, and finally applying a rigid body transformation to move this created cylinder, 
which is initially aligned on the z-axis, towards its desired 3-D location. For generating 
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a 3-D flexible tube called generalized cylinder, the method of the invention comprises 
starting from a straight cylinder, which is aligned on the z-axis, and which has a length s 
equal to the total length of the 3-D target path P. Then, the method comprises elastically 
warping this cylinder in order to fit the given 3-D path. 
5 Referring to FIG.l, the present method first comprises steps of: 

computing a 3-D path related to the tubular object of interest; for 
example the 3-D path is the centerline of a tubular object of interest to segment; 
defining segments on said 3-D path; 

creating an initial straight deformable cylindrical mesh model, of any 
10 kind of mesh, with a length defined along its longitudinal axis equal to the length of the 
3-D path; 

dividing this initial mesh model into segments of length related to the 
different segments of the 3-D path; and 

computing for each segment of the mesh, a rigid-body transformation 

15 that transforms the initial direction of the mesh into the direction of the related segment 
of the 3-D path, and applying this transformation to the vertices of the mesh 
corresponding to that segment. 

If the starting location of the straight cylinder is chosen to be at z-coordinate 0, 
then the z-coordinate of each cylinder point can be used as a correspondence with the 

20 chord-length parameterization of the 3-D path. For each location p(s) of the path, the 
direction u(s) of the local tangent is determined. Then, the cylinder points having a z- 
coordinate equal to s are warped into the plane issued from p(s) and orthogonal to u(s). 
These basic operations are illustrated by FIG.l. 
However, some artifacts might appear. 

25 First, if the 3-D path is not smooth, for example because large angles 

exist between two consecutive segments of the poly-line path, then the warped cylinder 
might cross itself, thus leading into an undesirable self intersecting mesh. Thus, a key 
problem is the potential apparition of self-intersections when the direction between two 
consecutive segments changes quickly. 

30 Then, when no care is taken with the choice of the rigid body 

transformation matching the z-axis into the u(s) orientation, this might lead to an 
undesirable torsion of the resulting mesh. Thus, another problem is the lack of 
continuity control, which can be seen as mesh torsion, during the transformation. 
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2) Limiting or avoiding self-intersections: 

A first cause of self-intersections is the rapid change of direction u(s) between 
two consecutive segments of the poly-line path. If the angle, denoted by 0, between the 
local tangent u(s) and the z-axis is plotted, a step-curve formed of steps as shown in 

5 FIG.2 is obtained. The steps are the horizontal segments A, B,.., represented in 

continuous lines in FIG.2. Since this angle G is directly used for computing the rigid- 
body transformation of each circle, as represented in FIG.l, reducing the effect of the 
steps helps preventing self-intersections. 

Referring to FIG.2, for reducing the effect of steps, the computed rotation for 

10 operating warping, corresponding to the segment "A" of the 3-D path, will be constant, 
and all the circles located between z=0 and z=l will be rotated with the same 3-D 
rotation. Similarly, all the circles corresponding to segment "B", i.e. having a z- 
coordinate in the range [1..2], will move with the same rigid body transformation. It is 
to be noted that, not only a 3-D rotation angle has to computed, but also a rotation axis, 

15 and similarly to the rotation angle, the rotation axis is piecewise constant along the s 
coordinate. 

2a) Limiting self-intersections by using linear blending: 

Self-intersections can be avoided if a unique transformation is not applied for 
each segment. Instead, the transformations are blended in between two consecutive 

20 segments using linear interpolation between two rotations. This linear interpolation 
from one 3-D rigid body transformation to another one can be easily solved using the 
formalism of Quaternions. Hence, according to the invention, in order to limit the self- 
intersections, the method comprises linearly blending the 3-D transformation between 
each segment, i.e. replacing the step-curve formed of segments A, B, etc, by the 

25 continuous curve represented by a doted line in FIG.2. This linear blending is achieved 
by considering the 3-D rotations as Quaternions (denoted by Qi, and Q2) with a 
quaternion of the form shown in (1): 

Q = cos - + (Lu x + j.u y + k.u z ) sin - (1) 

30 where 0 is the rotation angle, and ( u x , u y , u z ) the rotation axis. The blending itself is 
then straightforwardly a normalized linear combination of the two Quaternions Qi, Q2: 
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Non-linear blending could be investigated. But, for the example of tubular 
organ segmentation, the linear blending was proven stable enough. 

FIG.3 A and FIG.3B show the effect of the rotation blending on a 3-D path 
containing only 8 points, and also having quite large angle change from one segment to 
5 the other. In FIG.3A, it can be seen that, without 3-D rotation blending, the different 
circles intersect at the junction points, such as at points la, 2a, 3a, and the generated 
simplex mesh contains some self-intersections. In FIG.3B, it can be seen that the linear 
blending of the rotations helps the different circles to being deformed smoothly from 
one direction to the following one, resulting in a much more regular mesh, as shown at 
10 points lb, 2b, 3b. 

2b) Avoiding self-intersections by using automatic radius 

modulation: 

Linear blending of 3-d rigid transformation from one segment to the other does 
not always suffice to avoid self-intersections. Clearly, such self-intersections also 

15 depend on the relation between the local curvature of the 3-D path and the desired 
radius of the created mesh. If the latter is larger than the local radius of curvature, 
knowing that the radius of curvature is inversely proportional to the curvature, thus it is 
small when the curvature is high, then self-intersections occur. Thus, even if a smooth 
evolution of the rigid body transformation along with the s-coordinate is assured by the 

20 above-described operation of linear-blending, some self-intersection might still appear. 
The relation that exists between the radius, denoted by r, of the initial cylinder, the 
distance separating two consecutive circles, and the curvature, denoted by c, of the 3-D 
path, might influence the creation of such self-intersections. Trying to warp a cylinder 
with a large radius r on a very bent path will certainly lead to some serious problems. 

25 Hence, it is desirable to automatically reduce locally the diameter of the cylinder in 
highly curved zones. According to the invention, the mesh radius is adapted 
automatically, based on the curvature and sample distance of the points and the desired 
input radius. The method for tubular mesh creation comprises modulating the radius of 
the cylindrical mesh according to the local curvature. A shrinking factor combined with 

30 the 3-D rotation is calculated, thus leading to a similarity transform, provided that the 
curvature of the 3-D path can be computed at any point. Since the invention is related to 
organs, it is assumed that the provided polyline is smooth enough to use simple 


WO 2004/063988 


11 


PCT/IB2004/000063 


approximations. The curvature c(t) at curvilinear coordinate s=t of the 3-D path is 
computed as: 


c(t) = 


dp 

d 2 p 

II— II 3 

ds 

as 2 

1*1 


(3) 


where the first and second derivatives of the path are estimated with formulations (4) 
and (5): 


* ||pt-lPt[ + ||ptPt7l| 

* (||pt-2Pt-l| + ||Pt-lPt|| + ||ptPt+l|| + |Pt+lPt+2||) 

This is only valid when the length of each segment tends to zero and when they 

are quite homogeneous. To overcome this small regularity problem and to avoid large 

10 variation of the estimated curvature, a Gaussian filtering of the curvature values can be 

applied. Finally the computed curvature value for each point of the 3-D path is used to 

apply a shrinking factor to the different circles of the cylinder. This shrinking factor, 

denoted by k, depends on the radius of the initial cylinder r and the estimated radius of 

curvature, equal to 1/c, of the 3-D path: 

15 k= flifr.c(t)<=4 (6) 

|l/(r.c(t)) elsewhere v J 

Optionally, as with the rotation blending above, the shrinking factor can also be 
blended. The radius modulation technique can use a linear blending based on: 

(l-u)c(t 0 ) + uc(ti) (7) 
with u = x = (t - t 0 )/(ti - 1 0 ), (8) 
20 or this radius modulation technique can use a bi-cubic spline interpolation can be used, 
where u is replaced by v so that: 

v-3t 2 -2-c 3 . (9) 
FIG.5 shows an example of bi-cubic interpolation on computed curvatures. The 
curvilinear coordinates s relates to the segment indices. The bold lines represent some 
25 sample values of radius of curvatures 1/c, where c is the curvature, computed at 

locations s = 0 through s = 5. And the curve represented by a fine line shows the bi- 
cubic interpolation of these sparse data. The advantage of using a bi-cubic interpolation 
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here in comparison with a linear one is that the function has a continuous derivative (the 
tangents at the junction points are horizontal), hence the variation of the cylinder radius 
will vary slower. 

After all these precautions, if self-intersections still exists, then automatic mesh 
5 repairing, smoothing with internal force of the active contour algorithm might be 
applied, as described in the introduction part in relation with the transformations 
described in the prior art. 

FIG.6A and FIG.6B illustrate the benefit of automatic radius modulation. 
FIG.6A and FIG.6B show an example of mesh creation from a path in the colon. In 

10 FIG.6A, where the mesh model is created without using radius modulation, the colon 
mesh self-intersects in the bottom part 10a, where the 3-D path is very tortuous. Also, it 
very difficult to visualize some regions of the colon that are hidden by the bends of 
other regions. In FIG.6B, where the mesh model is created using the radius modulation 
algorithm, the self-intersections are largely reduced as shown in 10b. However, the 

15 general shape of the colon is not perturbed in the regions of restricted radii. In the other 
parts such as 1 la, the radius is unchanged between FIG.6A and FIG.6B. In the regions 
of restricted radii, the visualization and following of the different regions of the colon is 
greatly improved. 

3) Imposing minimal torsion: 

20 Mesh torsion is minimized when the distance between two consecutive 

- rotations, L e. rigid-body transformations, is minimal. The 3-D rotation is computed as 
the minimal rotation from the initial mesh direction, which is the z-axis, to the target 
segment u(s). This rotation is defined with an axis R = z A u(s), and with angle 9 = 
angle(z,u(s)). According to the invention, the parameters constituted by the axis R and 

25 the angle 9 are computed iteratively from one segment to the other. Thus, this method 
comprises computing the new rotation for a current segment u(s+l) as a composition of 
the found rotation for the previous segment u(s) and the minimal rotation from the 
previous segment u(s) to the current segment u(s+l). In fact, there is no unique 3-D 
rigid body transformation that maps the z-axis into the u(s) direction, which is the local 

30 tangent of the 3-D path. If a translation is needed, such a translation is easy to compute. 
But, the core problem consists in finding a 3-D rotation transforming the z vector 
(0,0,1) into u(s) = (u x/ u y , u z ) . The simplest rotation is the rotation around axis: 
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R = ZaU(s), (10) 

which reduces to R = (-u y , u x/ 0) , with a rotation angle 0 equals to: 

0 = -arcros(u z /|p(i)J). (11) 
This rotation, which may be called for simplicity: 
5 kl^s)) (12) 

maps the z vector to the u(s) direction, but any other rotation composed of a rotation 

around the z-axis (leaving the z-direction untouched) followed by R(z,u(s)) does it as 
well. Nevertheless, only one rotation fulfills the condition of minimal displacement and 
it is exactly R(z, u(s)) . The problem is that it must be taken care about the fact that no 
10 minimal torsion is to be found between z and u(s) , but rather between each step going 

from u(s) to u(s + 1) . Thus the rotation to be computed is not as easy as R(z, u(s)) , but 
corresponds to the incremental compositions: 

R(s) = R(i, m) o RM0), SO)) • ■ • ■ o R(u(s-1), ^)) , (1 3) 

where R(a,b) designates the minimal rotation mapping a into b . 

15 FIG.7A and FIG.7B illustrate the use of imposing minimal torsion. FIG.7A 

shows an example of mesh creation using only the minimal rotation between the z-axis 
and u(s). FIG.7B shows an example of mesh creation using the minimal rotation 
between the z-axis and u(s) and then the incremental rotation R(s) leading to a minimal 
torsion. In FIG.7A, it can be seen that torsion appears on the mesh because the cells are 

20 twisted around junction points, for example in regions 4a and 5a. Instead, in FIG.7B, the 
cells are kept aligned all over the mesh, ha particular, the cells are well aligned in 
regions 4b and 5b corresponding to the regions 4a and 5a of FIG.7A. 

On a tortuous synthetic example, it can also be observed that the minimal 
torsion algorithm gives fine results. FIG. 8 A represents an initial tubular mesh model 

25 from which a tortuous synthetic example is created and represented in FIG.8B and in 
FIG.8C under a different angle of view. FIG.8D represents the same example, with a 
finer resolution of the 3-D path. In FIG.8D, 200 segments were used in the 3-D path 
while in FIG.8B and FIG.8C only 50 segments were used in the 3-D path. In FIG.8A to 
FIG.8D, a straight line has been tagged in white on the initial cylinder 10, allowing to 

30 assessing the torsion on the final mesh in following the corresponding lines 1 1, 12, 13. 
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4) Desired initial conditions 

The above described method works with different kinds of 3-D paths. 
However, the best results are observed when no sharp angles are present. Hence, it is 
better to preliminary smooth the input 3-D path using any smoothing technique known 
5 to those skilled in the art. Still better results are also obtained when the segment lengths 
of the path are homogeneous. 

5) Medical viewing system and apparatus 

The above-described steps can be carried out by the viewing system of the 
invention. Fig.9 shows the basic components of an embodiment of an image viewing 

10 system in accordance to the present invention, incorporated in a medical examination 
apparatus. The medical examination apparatus 151 may include a bed 1 1 0 on which the 
patient lies or another element for localizing the patient relative to the imaging 
apparatus. The medical imaging apparatus 151 may be a CT scanner or other medical 
imaging apparatus such as x-rays or ultrasound apparatus. The image data produced by 

15 the apparatus 151 is fed to data processing means 153, such as a general-purpose 
computer, that carries out the steps of the method. The data processing means 153 is 
typically associated with a visualization device, such as a monitor 154, and an input 
device 155, such as a keyboard, or a mouse 156, pointing device, etc. operative by the 
user so that he can interact with the system. The data processing device 153 is 

20 programmed to implement a method of processing medical image data according to 
invention. In particular, the data processing device 153 has computing means and 
memory means to perform the steps of the method. A computer program product having 
pre-programmed instructions to carry out the method may also be implemented. 

The drawings and their description herein before illustrate rather than limit the 

25 invention. It will be evident that there are numerous alternatives that fall within the 
scope of the appended claims. Moreover, although the present invention has been 
described in terms of generating image data for display, the present invention is 
intended to cover substantially any form of visualization of the image data including, 
but not limited to, display on a display device, and printing. Any reference sign in a 

30 claim should not be construed as limiting the claim. 


